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Spatial manipulation of current flow in graphene could be achieved through the use of a tilted 
pn junction. We show through numerical simulation that a pseudo-Hall effect (i.e. non-equilibrium 
charge and current density accumulating along one of the sides of a graphene ribbon) can be observed 
under these conditions. The tilt angle and the pn transition length are two key parameters in tuning 
the strength of this effect. This phenomenon can be explained using classical trajectory via ray 
analysis, and is therefore relatively robust against disorder. Lastly, we propose and simulate a three 
terminal device that allows direct experimental access to the proposed effect. 

PACS numbers: 



I. INTRODUCTION 

A semiconductor pn junction where both sides of the 
junction are biased such that their Fermi surfaces are 
identical could potentially serve as an electron analogue 
of the Veselago lens pQ. Graphene, a zero band-gap two- 
dimensional semiconductor with Dirac-type linear energy 
dispersion [2j [3] [4] [5] , is an ideal medium for realizing this 
physical analogy. As recently advocated by Cheianov et 
al. [6], an abrupt and symmetrically biased graphene pn 
junction could function as an electron focusing device. 
An electronic superlattice of cascading pn junctions could 
serve as an electron beam colliminator as elaborated by 
Park et al. [7 j. In this letter, we propose utilizing a 
tilted pn junction to manipulate the current flow such 
that charge carriers preferably propagate along one edge 
of the sample using a set-up as illustrated in Fig. [T^i. This 
is achieved by controlling the following interface proper- 
ties: (i) tilt angle 5 and (ii) the extent of the pn transition 
region, D. The possibility to manipulate the spatial dis- 
tribution of the current density in graphene opens the 
door for novel electronic device concepts. 

Experimentally, graphene pn junctions are created 
through electrical means via a top/bottom gating scheme 
[lOj [TTJ [12]. Carrier transport across a conventional 
graphene pn junction exhibits highly angular selective 
behavior p~2j [13] [14j [15]. For example, in a symmetric 
pn junction, the transmission probability in the absence 
of magnetic field is given by T (fc m ) « e~ 7Tk ^ D/2k f [13], 
where fc// m is the Fermi and transverse wave- vector re- 
spectively. When km = 0, the transport across the pn 
junction would be reflect ionless, an hallmark of Klein 
tunneling [14]. Therefore, by geometrically tilting the 
graphene pn junction at an angle S as shown in Fig. [I] 
one expects that the maximum transmission now occurs 
for the transverse mode k m ~ kfsin5. Analogous to 
this physical situation is the problem of transport across 
a conventional pn junction in the presence of a mag- 
netic field, B. In the latter case, one uses the Lorentz 
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FIG. 1: (a) Schematic of a tilted pn junction device built on a 
graphene ribbon. The top and bottom gate allows the tunability of 
the electron/hole carrier density on each side of the junction. Tilt 
angle defined as 5. (b-e) Polar plots of the carrier transmission 
probability across a symmetric graphene pn junction for different 
values of magnetic field using the WKB model outlined in 8 (solid 
lines). The calculation is done for an experimentally typical pn 
junction with transition width of lOOnm and a built-in potential of 
OAeV i.e. ef = ±0.2eV for the n/p regions 9 . Similar plots for 
different 5 at B = OT using WKB are also shown (dashed lines). 



force to modify the carrier's trajectory. In the limit 
of large device width, one can impose the usual peri- 
odic boundary conditions and express the eigenstates as 
tf m (r) = e ik ^ y ip(x,B,k m ) [16 . The WKB transmission 
probability, Te(fc m ), in the presence of a B field is derived 
by Shytov and co-workers |5]. Fig. [l]>e plots Te(/u m ) 
for different B field strength. Reflectionless transmis- 
sion now occurs for the mode k m « kfsin0B, where 
Ob = sin~ 1 (vfB/E). The polar plots exhibits the char- 
acteristic leaf-shaped feature which rotates in the pres- 
ence of magnetic field. The thickness of the leaf defines 
the angular bandwidth (O) of the pn junction. Decreas- 
ing with increasing magnetic field is responsible for 
the degradation in conductance observed recently in ex- 
periments [9 . Suppose the device is large and the ef- 
fects from the boundaries are negligible, the transmission 
probability through a tilted junction could be described 
by a simple coordinate transformation, i.e. Tb(0 m — 5), 
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FIG. 2: NEGF derived mode-to-mode transmission probability function logi [T neg f (# m ,# n )] for a symmetric pn junction device biased 
at ef = OAeV, for tilt angles of 5 = 0°, 15°, 30°, 45° respectively. The device width is lOOnra and the pn transition length D = lOnra. 



as depicted in Fig. [l]>e. Both cases exhibit the sig- 
natures of a transverse current, i.e. Hall current in the 
magnetic field case. In this paper, we want to address 
the question, 'Could one engineer a pseudo-Hall effect in 
a graphene waveguide via a tilted pn junction as shown 
in Fig. [T^?'. 

II. THEORY AND METHODS 

The theoretical model we employ in this work is based 
on the Landauer-B'uttiker formalism where the transmis- 
sion function is computed within the framework of the 
non-equilibrium green function method [T7J [18] . The de- 
vice Hamiltonian is described within the tight binding 
formalism [H |20], 

H = YiiVia\di + Yiijta\dj (1) 

where a\/di are the creation/destruction operator at each 
atomic site i. Vi is the on-site potential energy, to be 
controlled by the top/bottom gates. Additional contri- 
butions due to local magnetization at the ribbon edges 
[2T] are not considered in this work, since our ribbon's 
width are relatively large. The open boundary condition 
for the quantum transport problem is embodied by con- 
tacts' self-energies, E s /^, solved using an iterative scheme 
outlined in [22]. The device Green function is then com- 
puted through, 

G(e f ) = (e f -H-Z s -Z d )- 1 (2) 

where e/ is the Fermi energy. However, direct matrix 
inversion of Eq. [2] usually proves to be computationally 
prohibitive. Therefore, one commonly resorts to recur- 
sive type techniques such as the recursive green function 
[23j[24] or the renormalization method [25] . In this work, 
we obtain the charge and current density of our device 
by combining familiar concepts from the recursive green 
function and the renormalization method. The detailed 
methodology is outlined in Appendix A. 

After solving for G(ef), we can compute the mode re- 
solved transmission probability function T™™^ at e/ via, 

T™ } = Tr [T s , m GT d , n G^] (3) 



where m/n denotes the modes in the source/drain con- 
tacts respectively. T s / d are known as the contacts' broad- 
ening functions which can be obtained from for each 
respective mode in the contacts i.e. T s / d = z2Jm(S s / d ). 
The mode-to-mode transmission function, T mn , is a use- 
ful quantity for analyzing the transport effects in the 
modal space in the presence of device non-homogeneities 
(see e.g. [26 ). Appendix B describes the procedure in 
obtaining the mode resolved contact self-energy S S)m in 
armchair edge graphene ribbon. 

III. RESULTS 

The graphene device that we investigate in this work is 
a lOOnm wide ribbon with armchair edges. For an arm- 
chair edge ribbon, the scattering states for an incoming 
source mode can be written as [27] . 

*™( r ) = 4ft x iMkm) -$K>(-km)}e ik * x (4) 
\J2W 

where (j)j^(k rn )=e' l1 ^' r e lkrnV and K/K' denotes the two in- 
equivalent Dirac points in the graphene's Brillouin zone. 
s(k) being the pseudo-spin, describing the A/B sublat- 
tice wavefunction. Since a graphene pn junction is anal- 
ogous to a negative refractive material in optics [6], it 
is useful to define an angular representation for the con- 
tact modes in order to facilitate discussion using sim- 
ple ray analysis. For the source/drain modes, we de- 
fine m / n = sin~ 1 (k m / n /kf) respectively, where n labels 
the modes in the drain. Fig. [2] plots the mode-to-mode 
transmission function for various tilted pn junc- 

tions (5 = 0°, 15°, 30°, 45°) biased symmetrically with a 
built-in potential of 0.8eV i.e. ef = OAeV on each side. 
When the pn interface tilt angle is zero, the solutions 
satisfy the 'Snell law' given by m = n (see Fig. ^f). 
When S ^ 0, the solutions are generally described by: 

a _ f m + 25 : K f . 

Un -\e m -2S : K' W 

as demonstrated in Fig. [2}d. Each source mode is an equal 
weight superposition of scattering states from K and K/ 
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FIG. 3: (a-c) Polar plots of the transmission probability 
T neg f(0m) = ZnT™ f for different 5, where T™ f is computed 
using the same device parameters as for Fig. [2] The WKB plots 
are obtained using the usual WKB formula |13]fsee text) 



valleys propagating in ±# m direction respectively. When 
£ 7^ 0, these two scattering states get scattered differ- 
ently, ending up in two different drain modes according 
to Eq. [5) However, when S exceeds a maximum tilt, i.e. 
5 maxi the solutions described by Eq. [5] fall outside of the 
available drain modes and new solutions emerge (see Fig. 
W). We find, 



\max (6 n ) 



(6) 



which can be easily deduced from Fig. [2] (for a lOOnm 
armchair edge ribbon, max(# n ) « 72°). We will revisit 
this point later. 

Summing over the drain modes, we obtained the trans- 
mission probability due to an incoming mode from the 
source: T negf (0 m ) = Z n T™ n gf . Fig. | is the polar plot 
of T neg f(6 rn ) for symmetric pn junction devices with dif- 
ferent 5. The WKB plots are obtained using the usual 
WKB formula [13] , and treating the two scattering states 
from each mode independently. When S increases gradu- 
ally from zero, the single leaf evolves into a doublet leaf 
structure as shown in Fig. ^$jp. This is because Klein tun- 
neling which originally occurs for the mode m ~ 0, now 
occurs for the modes m « ±5. For the latter, one would 
expect the maximum transmission probability to be ~ \ , 
since only one of the pair of scattering states from each 
mode satisfy the condition for Klein tunneling. How- 
ever, the NEGF result deviates from this simple picture, 
showing a notably higher maximum transmission than - . 
The reason for this discrepency is due to multiple scatter- 
ing with the sidewall i.e. sidewall enhanced transmission 
(SWET). With further increase in 5, this doublet leaf 
structure evolves into some triplet leaf feature (new solu- 
tions arise when S > 5 max ) and eventually the polar plot 
becomes noisy (not shown). Increasing of 5 extends the 
physical longitudinal distance of the tilted gate, thereby 
enhancing the mixing of the various transverse modes. 



A. Junction conductance 

Fig. |4^l plots the pn junction normalized conductance 
(a/a Q ) as a function of tilt angle 5 for different values of 
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FIG. 4: (a) Normalized conductance, cr/cro, of a pn junction as a 
function of the tilt angle 5, computed for various values of tran- 
sition length D and biased symmetrically at 6f = OAeV. ctq is 
defined as the conductance when 5 = 0, which by construction 
means cr/cro = 1 when 6 = 0. (b) Absolute conductance for the 
D = lOnm device with (i) perfectly smooth armchair edges (round 
symbol), roughened armchair edges i.e. RMS roughness of 1 atomic 
layer (square symbol) and perfectly smooth zigzag edges (triangle 
symbol). Dashed line is the result from simple WKB model. 



D, where <j is the conductance when 6 = 0°. The fol- 
lowing key observations can be made; (i) a/ao exhibits 
an initial increase with 5 and then decreases prominently 
when 5 exceeds a threshold angle, herein denoted as 5 t h, 
and (ii) the occurence of 5 t h can be delayed by employing 
a larger D. We also checked that the same trends hold 
true for devices with different widths, e/ and edge config- 
urations (i.e. zigzag edge ribbon shown in Fig. [4)3). The 
initial increase in conductance with S is a SWET phe- 
nomenon whose signature becomes more prominent with 
larger S. We shall discuss the plausible explanation for 
the existance of 5th, beyond which a /<jq degrades. In the 
transmission polar plot (Fig. [T]), increasing 6 rotates the 
leaf by a similar amount. The threshold of conductance 
degradation occurs at large enough S such that some of 
the states within the angular bandwidth would be back 
reflected into the source i.e. 



0th 



\ - n) 



(7) 



where Vt is the angular bandwidth of the transmission 
function. Since is larger for smaller D, S t h is also 
smaller. This explains the general trend we observe in 
our numerical calculation in Fig. [4j Let us compute Sth 
for the set of results in Fig. [4^i. Defining the Q to be 
the bandwidth where transmission probability is > 0.5, 
we have tt w 90°, 75°, 49°, 31° when D = 0,2.5, 5, lOnra 
respectively. This yields S t h ~ 45°, 53°, 66°, 74° respec- 
tively, in good agreement with the numerical result we 
obtained in Fig. [4^i. 

For device with D = lOnra, the junction conductance 
at S = 70° can exceed twice its value at S = 0°, as shown 
in Fig. [4^i. From a device perspective, this means that 
one could deliberately tilt the interface angle to enhance 
the on-state current. This could be useful for engineering 
a band-to-band tunneling transistor [28]. Next, we exam- 
ine the robustness of this effect in the presence of sidewall 
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FIG. 5: (a-c) shows the non-equilibrium real space longitudinal 
current density for pn devices with D = Onra, 5nm, lOnra respec- 
tively. All devices considered here have 5 = 30° and €f = OAeV. 
(d-e) are the same device as (c) excepts with sidewall roughness 
and pn interface roughness respectively. The RMS roughness are 
1 atomic layer and 1.7nm respectively, (f) is the device with 
D = lOnm and 5 = 45° with no disorder. 



disorder. Fig. plots the absolute junction conduc- 
tance of the D = lOnm device for the case with a per- 
fect sidewall and one where the sidewall exhibits a RMS 
roughness of one atomic layer. Evidently, the SWET phe- 
nomenon is highly sensitive to the characteristic of the 
sidewall. Therefore, chemically derived graphene ribbons 
[29] with smooth edges are needed to experimentally ob- 
serve these large conductance modulation with tilt angle. 
Device widths of the same order as the carrier's phase co- 
herence length is also required for SWET to occur. 
From Fabry Perot experiments [30], « lOOnra is ex- 
pected. 



B. Spatial current distribution 

Fig. [5^i-f shows the non-equilibrium spatial current in- 
tensity plots of a tilted graphene pn junction for different 
transition lengths and tilt angle. The current profile pop- 
ulates preferentially along the sides of the device, anal- 
ogous to the Hall effect. The following key observations 
can be made about Fig. [5^i-f; (i) edge populated current 
is observed to be more prominent with increasing D and 
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FIG. 6: (a) Schematic of a multiplexer device where the drain 
contact is partitioned into two via graphene stubs, (b) Ratio of 
the current through the two drain terminals as a function of S. 
Parameters assumed are D = lOnm, ef = OAeV and under small 
source drain bias. 



(ii) this effect can be enhanced by increasing 5 until 5 
exceeds a certain angle (see Fig. [5]p). The former effect 
is attributed to the suppression of normal modes current 
(i.e. T neg f(0 m « 0)) as a result of larger D. 

The appearance of edge populated current is a direct 
consequence of the negative refractive index property of 
the graphene pn junction. Each propagating mode from 
the source follows a classical trajectory into the drain 
contact according to Eq. [5] For example, when 5 = 15°, 
the incoming source modes m « ±15° would contribute 
the most current (see Fig. [3J3 and Fig. [2]o) . The mode 
m ~ 15° would end up in the drain modes 9 n ~ 45° (K) 
and n ~ — 15°(K / ), where the latter scattering state has 
a larger current contribution. (...) indicates the valley 
in which the incoming scattering state is residing. On 
the other hand, m ~ —15° would end up in the drain 
modes n « 15°(K) and n « -45°(K / ), where the for- 
mer has a larger contribution. Both scattering states, 
n w — 15°(K / ) and n « 15° (K), are propagating in the 
same direction, since K = — K 7 . This leads to the effect 
of pseudo-Hall current. However, when 8 > S max , new 
scattering states which do not follow the classical trajec- 
tories described by the 'Snell law' (i.e. Eq. [5| arises. 
They are responsible for overwhelming the edge popu- 
lated currents, resulting in their disappearance for the 
device with 6 = 45° (see Fig. |5f). Unlike the SWET, the 
phenomenon of edge populated current is fairly robust 
against various disorder such as edge roughness and pn 
interface roughness as shown in Fig. [5]i-e. One should 
be able to measure this effect experimentally [35] . 

We consider a possible experimental setup that allows 
direct access to the proposed effect as shown in Fig. [6^l. 
The drain is partitioned into two contacts through an 
upper/lower stub, with currents denoted as I1/I2 respec- 
tively. Fig. ^jp plots the ratio I1/I2 as a function of 5. 
The current asymmetry has an optimum value of 300% at 
S w 25° as shown. Through further device optimization, 
one should be able to engineer a device with a larger 
asymmetry ratio and achieve a semi-unipolar behavior 
through I2. 

In conclusion, we had peformed a numerical study 
of a tilted graphene pn junction, provides a detailed 
physical understanding of its transport properties, and 
highlighted the possibility of manipulating the current 
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to flow along the edges of the waveguide. 
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APPENDIX A: RENORMALIZATION AND 
RECURSIVE METHODS 
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This appendix documents the procedure we used for 
the computation of spatial charge/current density pro- 
files in the device. We consider a graphene ribbon with 
armchair edges as illustrated in Fig. [7] Device is infinite 
along x, the transport direction, and each supercell is 
represented by the dotted rectangular box. The Hamil- 
tonian, H, describing the graphene ribbon is formulated 
by treating only the nearest-neighbor interaction between 
the pz-orbitals p~9j [20] . Usually, this coupling energy is 
assumed to be t c = 3eV. By the same token, the su- 
percell would only interact with the adjacent supercells. 
The interation of a supercell with its neighboring cell on 
the right/left is represented by t/t^ respectively, while 
the intra-cell interation is denoted by a. r and a are 
matrices of size n s x n s , where n s is the number of basis 
functions in a supercell. H is the sum of these coupling 
energies and the electrostatic potential U(r). 

From a practical point of view, we are only interested 
in the scattering solutions within the central device do- 
main, denoted by Qq. Let Hl, Hq and H R be the Hamil- 
tonian description of ^l, and Qr respectively. The 
interaction between the Hl and Hq block is denoted by 
f . Through simple algebras, one can write the retarded 
green function at ef (Fermi energy) in fl G as follows [17] , 

G(r, r') = ((e f + irj)I -H -X s - = A' 1 (Al) 

£ s /d are known as the contact retarded self-energy and 
are defined as follows, 



= r ] g r L f 



9 r L = ((e f +if,)I- 
9r = (0/ +irj)I- 



Hr) 



(A2) 



The numerics for Eq. | A2| immediately becomes tractable 
when one notice that we only need the elements of g r L j R 
which are adjacent to the Q R D Qo an d &l H bound- 
aries. These 'surface elements', which are a smaller sub- 
set of g r L , R , are usually denoted by the surface green func- 
tion matrices g s L j R of size n s xn s . An iterative scheme is 
commonly used to compute g s L j R [22 . Once the contact 
retarded self-energies are determined, the device green 
function in Eq. |A1| can be computed by directly in- 



verting the matrix A, if computational resource is not 
a limiting factor. However, it usually prove to be compu- 
tationally prohibitive. Therefore, one commonly resorts 



FIG. 7: Schematic on a graphene ribbon with armchair edges. 
Each slice of supercell consist of an intra-cell interation denoted by 
a and a right/left neighboring cell interaction represented by r/r^ 
respectively. 



to recursive type techniques such as the recursive green 
function [23] [24] or the renormalization method [25] . In 
this appendix, we outlined a methodology which com- 
bines familiar concepts from the recursive green function 
and the decimation method in the computation of the 
various non-equilibrium green functions, from which we 
can obtain the charge and current density of our device. 

Suppose that we are only interested in the real space 
resolved charge and current density for the supercell j 
as shown in Fig. [7| The first step involves getting rid 
of the slices s = 2, 3, . . . , /i, /, . . . , n — 2, n — 1 from the 
system of equations stipulated in Eq. |A1[ m ade possible 
by the trigonal nature of matrix A. Eq. |Al| now becomes 
OA = I, where the LHS of this equation is explicitly 
written as, 



gn gu gij giu gin 

gn gu g^ g%k 

9jl 9ji 9 n 9jk 9jn 

gki gn gkj gkk gkn 

gnl gni gnj gnk gnn 



an an 





an 



uj 3% 








ajj ajk 

Q>kj dkk Q>kn 

&nk ttnn 



(A3) 



In this work, the matrix elements of A are systemati- 
cally derived. The elements not affected by the deci- 
mation process are ajj = [A] J -, aij = [A]p aji = [A]?, 



ajk = [A] J k and a^j = [A]j (the upper/lower index de- 
notes row/column respectively), an and an are obtained 
through a set of recursive formulae. We began with the 



initialization a% = [A] v 
mulae for an and an are, 



and a n = 



r^. The recursive for- 



pi 



u 



([Ay-izl-rpr^y 1 



a il Pi T 



(A4) 



where p] 1 - 



ai- 3 and 



0. The desired solutions are a? 
an = a J n ~- Note that a different set of recursive formula 
is needed if the intercell coupling is different for each 
supercell. akk and a^ n are obtained through a similar set 
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of recursive formulae. We began with the initialization 
a kk = [A] k and a kn = r. The recursive formulae for akk 
and dkn are 7 
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a kk 
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= ([At-l-l-Spr^y 1 (A5) 


9jj = 


where pf 
and a kn 


= 0. 

n— 
= a kn 


The desired solutions are akk = a kk^~ 2 
3 ~ . an and an are obtained through 


9ji = 
9ki 



a similar set of recursive formulae. We began with the 



initialization ai 
formulae are; 



o _ 
i 



[A]{ and aj i = r. The recursive 
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where qf 
and aii = 



u 



([Ar 



(A6) 



The desired solutions are an = a J u 3 
a nn and a n /c are obtained through a 
similar set of recursive formulae. We began with the 
initialization a^ n = [A]™ and aP nk = . The recursive 
formulae are; 



a nn 

u 
a nk 

q u r 



a nn a nk °r \ a n k ) 
u-l u f 
a nk Qr T 



n — u 
n — u 



rq r 



(A7) 



a n ~ 3 ~ 



where = 0. The desired solutions are a nn 
and a n k =Jh ] h J ~ • Performing the recursive procedure 
in Eq. A4||A7[ one can then obtain the full information 
of the matrix A. 

We are now ready to compute the charge and current 
density for the supercell j. The key quantity is the elec- 
tron correlation function given by, 



where S m = + are known as the in-scattering 
self energies. It is given by, 



^s/d 



;/ s/d (s s/d -£t /d ) (A8 ) 



where f s / d are the Fermi occupation factor in the source 
and drain contacts. 

The key concepts in recursive solution of G can now be 
employed to solve Eq. |A3| Specifically, we only require 
the solution to gjj, gj n , gk n , 9ji and gki for reasons that 
would be apparent later. We would need the following 
recursive formulae, 



where Q° = and it yields us ^ 5 = g n 
arrive at the following results, 



-tf[A]lg nn 
-n 3 [A]lg kn 

tf-n 4 [AM n 

n 3 - n 3 [A}l(n 3 [A]lg kk ) 

(-i) 2 (^]^ 2 »,,) 

2o3r 



We can then 



= (-1) 3 (0 1 [A]^0 2 [l]20 3 [A]^ fe ) J (A10) 



With these block elements information of G, we are now 
ready to compute the charge and current density. 

We are interested in the electron density, n(r), of the 
supercell j given by, 



[G 



n]J 
j 



[GK[sr];[Gt]; + [G]i[s™]:[Gt]; 
5il [Er]; 5 ] 1+5 ,„[E»]>]„ (ah) 



where we had make use of the fact that ^ l J^ d are non- 
zero only for j = l,n slices respectively. The current 
density, j(r), flowing between the supercell j and j + 1 
is computed via, 



j(r) 



2q 
h 



([^+i [G n \i + 1 - [A]] + 1 \G n Y j+l ) (A12) 



where, 



[G»g +1 = [G}{ +1 [ZT]\[G^) + [G}i +1 [Xf]l[G^ 



= g kl [Ht]\g),+9kn [V?]y ir , 



(A13) 



[G}i[K n ]lm] +1 +[Gi[^x[ G T j+1 

fti^lIflL+^^Hi (A14) 



Therefore, we have completed our procedure in comput- 
ing the charge and current density for the supercell j. At 
any time, our numerical procedure only requires direct 
matrix inversion of size n s x n s , where n s is the num- 
ber of basis functions in a supercell. By parallelizing the 
computations, we can compute the charge and current 
density for any number of supercells within the device 
domain. In our work, we had computed for a dozens of 
supercell to give us the required spatial resolution of the 
charge and current density related graphical plots in the 
main paper. 



APPENDIX B: MODE RESOLVED CONTACT 
SELF ENERGY 



[G\l = W-W[A)l 



[Gf 



q+r 



-W[Af q+1 [Gf +1 



q+r 
v+1 



(A9) 



The mode-to-mode transmission function, T 71 



is a 



useful quantity for analyzing the transport effects in the 
modal space in the presence of device non-homogeneities 
(see e.g. [26 ). It is given by T mn = Tr [T s ^GT d ^] , 



7 



where m/n denotes the modes in the source/drain con- 
tacts respectively. T s / d are known as the contacts' broad- 
ening functions which can be obtained from S s /^ for each 
respective modes in the contacts i.e. T s / d = i2Im(T is / d ). 
This appendix describe the procedure in obtaining the 
mode resolved contact self-energy H s , m in armchair edge 
graphene ribbon. 

In armchair ribbon, the analytical solutions of the 
wavefunction and energy dispersion is known analyti- 
cally [51] . One could construct a unitary operator V 
which perform the transformation from real space to 
mode space. Zhao and Guo [32] outlined the recipe for 
doing so. For mode m, its propagation along the lat- 
tice chain could be described by an on-site and coupling 
matrix a and /3 respectively, 



a = 



U m 

t c 

t c LU m 

Um 





- 








o- 


p = 




























Ac 








0_ 



(Bl) 



of carbon layers along the width direction. In the pa- 
per, the angular representation for mode m is given by 
m — sin~ x [{t c — Lu m )/ef]. The self-energy for the semi- 
infinite leads of this lattice chain is denoted by S m and 
could be computed rather inexpensively via the usual 
technique described in [22 or analytically as discussed in 
[32] . Finally, the real space form of the mode-resolved 
self-energy is given by, 



v (3 ro ® i m ) y f 



(B2) 



V is a n s x n s unitary matrix, whose elements values 
are assigned as described in [32 . I m is a n s /4 x n s /4 
matrix with elements given by I m (i,j) = <Si >m 5j jm . £ s , m 
is therefore of size n s xn s . To ensure that the procedure 
is correct, we check the following sum rule, 



S s — E s i + E. 



3.2 T ^s,3 



(B3) 



where uo m = 2t c cos (ra7r/(2L + 1)), L being the number This completes the objective of this appendix. 
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